
x = 1:5; %  idem: x  = linspace(1,5,5)
y = 1./x;

z = [];
for j = 1:101
	z = [z,24/25 + j/25];
endfor

w = lagrange(x,y,z);

plot(z,1./z,'b',z,w,'r');
hold off




